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We consider a suspended elastic rod under longitudinal compression. The compression can be 
used to adjust potential energy for transverse displacements from harmonic to double well regime. 
The two minima in potential energy curve describe two possible buckled states. Using transition 
state theory (TST) we have calculated the rate of conversion from one state to other. If the strain 
e = 4e c the simple TST rate diverges. We suggest a method to correct this divergence for quantum 
calculations. We also find that zero point energy contributions can be quite large so that single 
Cn ' mode calculations can lead to large errors in the rate. 
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INTRODUCTION 



Considerable attention has recently been paid to two-state nano-mechanical systems [1|-|7| and the possibility of 

i 1 observing quantum effects in them. Roukes et al. [l| proposed to use an electrostatically flexed cantilever to explore 

the possibility of macroscopic quantum tunnelling in a nano-mechanical system. Carr et al. [5j, [8J suggested using the 

two buckled states of a nanorod and investigated the possibility of observing quantum effects. A suspended elastic rod 

of rectangular cross section under longitudinal compression is considered. As the compressional strain is increased to 

qj • the buckling instability [8], the frequency of the fundamental vibrational mode drops continuously to zero. Beyond 

_^ the instability, the system has a double well potential for the transverse motion (sec Fig. Q]). The two minima in the 

potential energy curve describe the two possible buckled states at that particular strain [8( and the system can change 

Y\ from one to the other by thermal fluctuations or quantum tunneling. Since both the well depth and assymetry are 

. — • tunable, a variety of quantum phenomena can be explored, including zero-point fluctuations, tunneling and coherent 

superposition of macroscopically distinct states. In this paper we have calculated the rate of transition between the 

pH ■ two thermally equilibrium states of a buckled nanorod/nanotube at fixed strain using quantum mechanical version 

Q_/ of transition state theory. Although this mechanical system has analogies to the superconducting interference device 

- ' ■ in which the first observation of a coherent superposition of macroscopically distinct states was reported |9j , but our 

method strictly address the thermally induced incoherent relaxation rate from one buckled state to the other at a 

£> ■ fixed strain. It may be possible by compressing a silicon nanorod or a carbon nanotube at low temperature to cause 

a crossover between the quantum and thermal fluctuation regime. In an earlier paper we have analyzed the classical 

*/"") ■ version of this problem 

en 

en 

(^ ■ THE MODEL 

o. 

We use L, w and d (satisfying L >> w >> d) to denote the length, width and thickness of the rod jg, Illl4l3| . 



F is the linear modulus (energy per unit length) of the rod and is related to the elastic modulus Q of the material 
by F — Qwd. The bending moment k is given by k 2 = d 2 /12 [8|, [l4| for a bar of rectangular cross section. We 



K^ take the length of the uncompressed rod to be Lq. We apply compression on the two ends, reducing the separation 
between the two to L. If y(x) denotes the displacement of the rod in the 'd' direction, then the length of the rod 
Ltotai = Jo dx-\/l + (y') 2 ~ L + 1/2 L dx(y') 2 . The compression causes a contribution to the potential energy 

Velastic = F / ' {2L ){L to tal — L Q ) . (1) 

In addition, bending of the rod in the 'd' direction cause bending energy 

.2 /r, / J„/„//\2 



V b = Fk 2 /2 / dx(y"y. (2) 

Jo 



Thus the total potential energy is given by 

V[y{x)] - 1/2 / dx{Fn 2 {y") 2 + Fe{y') 2 ) + F/(8L )( f dx{y') 2 ) 2 + F/(2L )(L - L ) 2 . (3) 

Jo Jo 

Here e = (L — Lq)/Lq is the strain. 
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FIG. 1: Potential energy V as a function of the fundamental mode displacement Y. The shape of the potential energy is 
harmonic for e > e c , quartic for e — e c = critical strain (e c < 0) and a double well for e < e c . 
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FIG. 2: The rod under compression: The central figure (A) shows the uncompressed rod of length Lq. On compressing to 
length L, the rod buckle, either to B- or to B+. 



EXTREMA OF THE FUNCTIONAL V[Y(X)] 
Extremisation of potential energy functional with respect to y(x) leads to 



(4) 



and the hinged end points have boundary conditions y(0) = y(L) — and y"(0) = y"{L) — (hinged boundary 
conditions are chosen for computational simplicity). If e > e c = —k 2 tt 2 /L 2 then, the only solution to Eq. Q is 
y(x) — if e > e c . But if e < e c , two bucked states are possible. They are 



y{x) = ±A^2/L sia(wx/L) , 



(5) 



with A = ^2LoL 2 (e c — s)/tt 2 . For e < e c , all the normal modes of vibration about these are stable. The solution 
y(x) — is now a saddle point. 



THE DYNAMICS 

One can calculate the barrier height for the process of going from one buckled state to the other over the saddle 
(linear geometry) as 

pLinear rpr /_ _ \2 



^m^Zr = FL (s s c ) 2 /2. (6) 



-L 



The kinetic energy of the rod is [i/2f Q y 2 dx, where [i — tti/Lq is the mass per unit length. Using the boundary 
conditions for the hinged end points, we find the normal modes of the rod, y(x,t) — y n (x)e zul " t . At the saddle point, 
we obtain 

y n {x) = A n y/2/Lsm(nivx/L), (7) 

with n = 1, 2, 3.... The normal mode frequencies at the saddle point are given by 

W Lnear,n = W n y/fl 2 - £ /e c , (8) 



where ujq — tt 2 k/L 2 ^F/ /i. n — 1 is the unstable mode and it has the imaginary frequency 



U 'Linear, 1 — z "Lineori [") 



where ^Linear — <^o V ' zjs-c — 1- For the buckled state, the normal modes are the same as at the saddle point, but 
the normal mode frequencies are different. They are 

Ui n = UoTiyn 2 — 1 torn > 1, (10) 



while uj\ = Wov2 {e/s c — !)• The rate expression using classical TST is (Eq. 3.14 of reference [l5() is 



JV 



n=l 



^classical q FT , , /, ,t „ — A-E„'""*" /fcT /ii\ 

H f - il ^near ll\U n /U Linearn e Ba '" e " , (.UJ 



where N denotes the total number of transverse modes of the rod. One makes a negligible error by taking the value 
of N to be infinity and this leads to [10( 



R f assical - 2^ f [ 2 - V z) r I 2 + V z ) (£c - £)e ~^ (e ~ r (12) 



The rate expression using quantum TST is given by 



^quantum _ ^Ltnear^ 2kT ^ 6c ' . u/ ^1 s, 77 Sm ^(^f) /, „\ 

«/ - ssc»s=) smh( »? ) M 8i „ ll( i%^j- (13) 



The above rates have a problem. As \Jeje c — > 2, uJ Linear 2 — > the second buckling instability sets in and the rates 
diverge. Also the quantum rate expression diverge at temperature T — h£l Linear / \2k-n) . 

BEYOND THE SECOND BUCKLING INSTABILITY 

As the rod is compressed, first the mode A\ \/2jL sin(nx/L) becomes unstable and this is the first buckling in- 
stability and the rod buckles as a result of this. The length at which this occurs shall be denoted by L/. If one 
supposes that the rod is compressed further keeping the straight rod configuration, then at a length L Sl the mode 
A-2. \J2JL sin(27ra/L) too would become unstable and this is the second buckling instability. For e > 4e c , there is 
only one saddle point. But for e < Ae c , due to the second instability, the saddle point bifurcates into two. In a similar 
fashion one can have the third instability at a length L t etc. In order to analyze the rate near and beyond the second 
buckling instability, we assume that the displacement has the form 

y (x) = A lx /2/Lsm(Trx/L) + A 2 \/2/Lahx(2nx/L). (14) 
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FIG. 3: Plot of AE^Z e r Zr against (L - L ) for a silicon rod of dimensions L = 1000 A, w = 200 A and d = 100 A. For this 
rod the first two buckling instabilities occur at Lf — Lq — —8.364 A and L s — La = —35.354 A respectively. The first buckling 
instability is shown by an arrow and the second buckling instability is not shown in figure. 
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FIG. 4: Plot of barrier height against (L — Lo), for a silicon rod of dimensions Lq — 500 A, w = 20 A, d— 10 A. For this rod the 
first three buckling instabilities occur at Lf —Lq — —0.1646 A, L s —Lq — — 0.6597 A and L t — Lo — —1.4893 A respectively. The 
first two buckling instabilities are shown in figure using arrows and the third buckling instability is not shown in figure. Solid 
line is for linear saddle point (valid in the regime L > L s ), dashed line is for bent saddle point (valid in the regime L < L s ). 



Using this, the elastic potential energy is given by 

Ftt 4 (Ai 



V(A 1 ,A 2 



4A|) 2 



8L 4 L 



FTr 2 Ai(e-e c ) 2Ftt 2 A 2 2 (e - 4e 



2L 2 



L 2 



(15) 



Finding the extrema of this potential leads to the following three solutions for (A\, A%): (a) (0, 0): this is the 
straight rod configuration. Between first and second buckling (i.e. L s < L < Lf ), this is the saddle point. But 
after the second buckling, it is no longer a saddle, but it becomes a hill top. It has the energy Ehuitop = 0. (b) 
(±2/7ry / LL (e c — e), 0): These are the buckled states and both of them have the same energy E^ = —FL (s — e c ) 2 /2. 
(c) (0, ±l/7T-y/LLo(4£ c — e)): These are the two new saddle points that arise from the bifurcation of the one that 




-1.60 -1.58 
(L-L J/Angstrom 



FIG. 5: Plot of AJSg^rrier against (L — Lo), for a multiwalled carbon nanotube of dimensions Lo = 5000 A, di — 50 A, 
d,2 — 100 A. For this nanotube the first two instabilities occur at Lf — Lo — —1.543 A and L s — Lo — —6.184 A rexpectively. 
The first buckling instability is shown by an arrow and the second buckling instability is not shown in figure. 
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FIG. 6: Plot of A_Ega"rier against (L — Lo), for a multiwalled carbon nanotube of dimensions Lq — 500 A, d\ — 10 A, cfe = 50 A. 
For this nanotube the first two buckling instabilities occur at Lf — Lo — —3.2497 A and L s — Lo = —13.5555 A respectively. 
The first buckling instability is shown by an arrow and the second buckling instability is not shown in figure. 



existed for 4e c < e. At these saddle points, the rod has a bent (S shaped) geometry. These two have the same energy 



E 



Bent 

■nldlc. 



-FL (e~4e c ) 2 /2. 



Beyond the second buckling instability, the barrier height is given by 



A£f^ ier = -3FL s c (-2e + 5e c )/2. 



Near the saddle, the normal mode frequencies are given by: 



^Bent.l - l ^Bent, 



(16) 
(17) 
(18) 



6 
wL.t.2=<W8(e/eo-4), (19) 



and for n > 2, 



4ent,« = wonV^ 2 - 4. (20) 

In the above ur Bent x has an imaginary frequency with ft Bent = v3^o- 

RATE BEYOND THE SECOND BUCKLING INSTABLILTY 

Now the classical rate beyond the second buckling instability can be calculated taking the saddle to be the bent 
configuration |l0| : 



R clas Sl cal = 4^/^ /i£__L e ^8? a (-2 e +fa e ) # (21) 

V 4e c - e 

It is interesting that the normal modes for this saddle retain their stability, irrespective of what the compression is. 
The quantum rate expression is given by 

quantum = ^Bentg ^ ^ sin h(-^) TT S|*£J , (22) 

™(^) 2fcT Msinh(%^) 

The above rates have a problem. As \/e/e c —¥ 2, w]j ent 2 — » the second buckling instability sets in and the rates 
diverge. Also the quantum rate expression diverge at temperature T — hQ Bent/ (2kir). 

RATE NEAR THE SECOND BUCKLING INSTABLILTY 



Near the second buckling instability {y/e/e c — > 2) vanishes, causing the rates in Eq. (fl"2j) . Eq. (|T3j) . Eq. (|21|) and 
Eq. (|22[) to diverge. The cure for the divergence is simple for the classical rate and is given below in the equations 
(j23|) and (|24[) . All that one has to do is to include the quartic term in A 2 of Eq. (|15[) in the evaluation of partition 
function for the second mode at the saddle. All the other modes (at the saddle as well as at the reactant) are treated 
as harmonic. In the regime where L > L s the rate is then given by (saddle is the straight rod) [10| 



FL„ 



^^^1 = £ — /<nt(4 _ L,-JL—) W^-p(f - i)r(3- ./f )r(3 + ,/f ), (23) 

where fi n t(a, b) = J_ dye . In the regime where L < L s the rate is given by (saddle is the bent rod) [10( 



FL 



£(e- E c) 2F{-s c ).e [3 , fA e hT 

L y /x e c V 7r e c 2FL Q ei 



Ri al = - r \/— V"^(- " l)\/zAnt(4- - ^rm)- ( 24 ) 



Now we follow the work of Voth, Chandler and Miller [16l | for rate calculation using quantum mechanical transition 
state theory near the second buckling instability (where 0J Linear 2 is small) . 

The QTST rate under the harmonic approximation diverges at y/e/e c — > 2, but using this method it is possible 
to avoid the divergence |16| . For this we go beyond the harmonic approximations for the first two modes. The 
Hamiltonian for the first two modes (A±, A2) at the saddle may be written as: 

H = p 2 l /(2n)+p 2 2 /(2n) + V(A 1 ,A 2 ), (25) 

where pi = fiAi is the momentum operator canonically conjugate to A%. The Lagrangian is 

L(A 1 ,A 2 ) = M /2(Ai 2 + A2) - V(A 1 ,A 2 ). (26) 
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FIG. 7: Plot of rate of crossing from one buckled state to the other using classical transition state theory, for a silicon rod 
of dimensions Lo = 1000 A, w = 200 A, d = 100 A at T — 300 K. For this rod the first two buckling instabilities occur at 
Lf — Lo = —8.364 A and L s — Lo — —35.354 A respectively. The first buckling instability is shown by an arrow and the second 
buckling instability is not shown in figure. 



The "centroid" partition function for these two coupled modes is given by the path integral 

Q*= f DA 2 (t) /"£Ui(T)<5(Ai-4)e~ S[Al(T),A2(T)1//i . 
Here Ai denotes the position of the centroid, defined by 



A, 



0h 



1/OSft) / drA^r). 



(27) 



(28) 



Note that the centroid for the first mode is constrained at A\, where A\ is the value at the saddle [16| and is equal 
to zero. We calculate the frequency of the unstable mode O and the partition function for the second mode Q\ using 
a variational principle based on the trial action la, Il7| 

dr{l/2 Ai A 1 2 (r) - 1/2^ 2 A\{t) + L 2 (A 2 ) + l/2^A 2 (r) + 1/2^(A 2 )(A 2 (t) - A 2 ) 2 }. (29) 



Atrial 



One can determine fi 2 , uj 2 {A 2 ) and L 2 (A 2 ) variationally, so as to get the best possible value for Q* . Once these values 
are obtained, we can proceed to calculate the rate, because in the trial action, the two modes are decoupled. Now the 
partition function of the second vibrational mode at the saddle Qjj; may be approximated by 



Ql = J(fj,kT)/2nh z 



W 2 (A 2 ) 

dA 2 e kT 



with an effective potential 



W 2 (A 2 ) = kT\og[2kT/(huj 2 (A 2 ))smh(hLj 2 (A 2 )/(2kT))}+L 2 (A 2 ). 



(30) 



(31) 



The tunneling current for the first mode may be taken as £1 / (2ir sm(i^7p)) , where fi is the variationally determined 
frequency of the reactive mode. In the regime where L > L s the rate may be calculated using (transition state is 
assumed to be straight rod) 
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(32) 



• 1 f T irirnr n \ 

n=3 smh( 2fcT ' ) 
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FIG. 8: Plot of rate of crossing from one buckled state to the other using classical TST, for a silicon rod of dimensions 
Lo = 500 A,w = 20 A, d = 10 A at T = 300 K. For this rod the first three buckling instabilities occur at Lf - L = -0.1646 A, 
L a — Lo — —0.6597 A and L t — Lq = —1.4893 A respectively. The first two buckling instabilities are shown by arrows and the 
third buckling instability is not shown in figure. Dotted line is for linear saddle point, dashed line is for bent saddle point and 
solid line include quartic terms for the second mode at saddle point. 



In the regime where L < L s the rate is given by (transition state is assumed to be bent rod) 

2kT' 



R 



quantum 



^-^'^\ in ^ )shM ^ )Qi fi^m 



Trsin(^) -v 2fcT'--- v 2fcT 

Also the quantum rate expression diverge at temperature T — hfl/(2kir 



(33) 



=3 sinh(- 



2feT 



RESULTS 



In the following we discuss the results of our calculation for rectangular silocon rods and cylindrical multiwalled 
carbon nanotubes. We first consider the Silicon rods of different dimensions, as summarized in Table. HI Si has an 
Young's modulus Q = 130 GPa and density p — 2230 kg.m~ 3 . First we consider a rod of dimensions 100 nm x 20 nm x 
10 nm, considered by Carr et al [8]. In Fig. [3] we plot the activation energy against the compression. At L = Lf, 
the activation energy is zero and as one compresses the rod, it increases rapidly, as it is proportional to (e — e c ) . 
Over a very short range of compression (~ 0.14 A) it increases to about 6 kcal/mol. The classical rate at T=300 K 
is plotted in Fig. At L = Lf, the potential is very flat and near this L the pre-factor (i.e., the attempt frequency) 
in the rate expression vanishes. This is the reason why the rate goes to zero as L — > Lf. It is found that the rate 
increases at first as one compresses. This is due to the increase in the pre-factor for the rate from zero. Then the rate 
decreases due to the increase in the barrier height. For the quantum calculation, we choose the value of N to be equal 
to the transverse degrees of freedom that would be there if one considered an atomistic model for the rod. Thus, for 
this rod, we took N = 425. We also report (see Table |T| the quantum enhancement factor, defined by r=Quantum 
Rate/Classical Rate. 

The values of T was found to be ~ 1 implying that quantum effects are not important for this rod, within the 
observable range of compression. For a rod of dimensions 2000 nm x 20 nm x 10 nm, plots of the rates at 0.01 K is 
given in Fig. [IT] and h shows that quantum effects lead to an increase in the rate by about a factor of 10. However, 
the overall rate rapidly decreases to very low values as one compresses the rod by about 0.01 A (see the Fig. [TT]) 
and hence it would be very difficult to observe this quantum enhancement experimentally. We have also performed 
calculations for a hypothetical rod of dimensions 50 nm x 2 nm x 1 nm. It is yet not possible yet to have Si rod of 
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FIG. 9: Plot of rate of crossing from one buckled state to the other using classical TST, for a multiwalled carbon nanotube of 
dimensions Lo = 5000 A, d\ — 50 A, d,2 = 100 A. For this nanotube the first two instabilities occur at Lf — Lo — —1.543 A and 
L s — Lo = —6.184 A rexpectively. The first buckling instability is shown by an arrow and the second buckling instability is not 
shown in figure. 
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FIG. 10: Plot of rate of crossing from one buckled state to the other using classical TST, for a multiwalled carbon nanotube of 
dimensions Lo — 500 A, ct = 10 A, di = 50 A. For this nanotube the first two buckling instabilities occur at Lf — Lo = — 3. 2497 A 
and L s — Lo — —13.5555 A respectively. The first buckling instability is shown by an arrow and the second buckling instability 
is not shown in figure. 



these dimensions. However it should be possible to synthesize molecular rods of these dimensions [18]. In Fig. [4] we 
have plotted barrier height as a function of compression. In the regime L > L s , the transition state has straight rod 
configuration and in the regime L < L s , the transition state has bent configuration. Fig. [8] shows the classical rate 
against compression, made at a temperature of 300 K. The rate obtained using the linear transition state is seen to 
diverge at the second buckling instability, but is finite for all L > L s . Similarly, beyond the second buckling instability 
the rate is calculated using the bent saddle. Close to the instability, this rate too diverges, but a well behaved rate 
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FIG. 11: Plot of logarithm of rate of crossing from one buckled state to the other using both classical TST and quantum TST, 
for a silicon rod of dimensions Lo = 20000 A, w — 200 A and d — 100 A at T — 0.01 K. For this rod the first two buckling 
instabilities occur at Lf — Lo = —0.41125 A and L s — Lo = —1.6452 A respectively. The first buckling instability is shown by 
an arrow and the second buckling instability is not shown in figure. The solid line is the result using classical TST, dashed line 
is the result using quantum TST. 
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FIG. 12: Plot of logarithm of rate of crossing from one buckled state to the other using both classical TST and quantum TST, 
for a silicon rod of dimensions Lo = 500 A,w — 20 A, d — 10 A at T — 0.01 K. For this rod the first three buckling instabilities 
occur at Lf — Lo = —0.1646 A, L s — L = —0.6597 A and L t — Lo = —1.4893 A respectively. The first two buckling instabilities 
are shown by arrows and the third buckling instability is not shown in figure. The solid line is the result using classical TST, 
dashed line is the result using quantum TST. 



can be calculated using the approach of Voth et al. 16] outlined above. In the quantum rate calculation we have 
taken contributions from 213 normal modes of this rod. We have compared the quantum rate with the classical rate 
at 0.01 K (Fig. [T2"]) . There is a quantum enhancement in the rate of roughly 10 6 . This occurs as one compresses the 
rod by about 0.015 A. Again, fabricating such a rod and doing an experiment under such conditions is a formidable 
challenge. 

Now we consider the case of carbon nanotubes used by Carr, Lawrence and Wybourne [8], as summerized in Table 
1. The larger tube have dimensions typical of the multiwalled nanotube whose vibrational properties were studied 
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FIG. 13: Plot of logarithm of rate of crossing from one buckled state to the other using both classical TST and quantum 
TST, for a multiwalled carbon nanotube of dimensions Lo = 5000 A, d\ — 50 A, d,2 — 100 A. For this nanotube the first two 
instabilities occur at Lf — Lq — —1.543 A and L 3 — Lo = —6.184 A rexpectively. The first buckling instability is shown by an 
arrow and the second buckling instability is not shown in figure. The solid line is the result using classical TST, dotted line is 
the result using quantum TST. The first buckling instability is shown by an arrow. 



by Treacy et al. [19j . The smaller nanotube is most probably the smallest nanotube that would support buckling 
and retain its elastic integrity [20J. In case of cylindrical tube, the bending moment k is given by k 2 = {d\ + d|)/16 
|8(, where d\ and di are inner and outer diameters. The Young's modulus and density of carbon nanotubes are [191 ] 
Q = 1.8 TPa and p = 215(%/m 3 . 

For a multiwalled carbon nanotube of dimensions Lq = 5000 A, d\ = 50 A, e?2 = 100 A. [8], plots of the rates at 
0.01 K is given in Fig. [13] and it shows that quantum effects lead to an increase in the rate by about a factor of 10. 
However, the overall rate rapidly decreases to very low values as one compresses the rod by about 0.001 A and hence 
it would not be very easy to observe this quantum enhancement experimentally. 

Now we consider a multiwalled carbon nanotube of dimensions Lq = 500 A, d\ — 10 A, di — 50 A [8(. In Fig. [6] we 
plot the activation energy against the compressions of the nanotube, it increases very rapidly. Over a very short range 
of compression (rj 0.1 A it increases to about 5 kcal/mol. The classical rate at T=300 K is plotted in Fig. [TUl The 
value of r is found to be rj 10 at T=0.01, implying that quantum effects are important at this temperature, within 
the observable range of compression (Fig. |T4")) . 



CONCLUSIONS 



We now ask, what is the origin of a high quantum enhancement factor in some of the cases? This is not due to 
tunneling of the reactive mode l^_ 8j], but is a zero point energy effect. The earlier analysis took only the reactive 
mode in their calculations }8l llll4l3j . but our analysis takes all the modes. This leads surprisingly to a decrease in the 
effective activation energy. The effective quantum mechanical barrier height may be written as AEg™^^ + AE zero , 
where AE zero is the difference in zero point energy of the modes between the buckled state and the transition state. 

For the linear transition state AE zero = Yl n =2 ^ lu, Q n ( \/ n2 ~ £ / £ c — Vn 2 — 1 ) = Yl n =2 ^ w o(-^ — 1) ( £ / £ c ~ 1) /2. 
Notice that AE zero is negative and hence leads to a lowering of the barrier height, which is proportional to N. This 
is the reason why the quantum effect is more pronounced for the 2000 nm bar (provided w and d are the same) as 
may be seen in Table |U The only possibility for observing quantum effects seems to be exciting the rod to a higher 
level in the buckled potential, suggested by Blencowe [7[. Thus, for the 2000 nm bar if one keeps L — Lo = — 0.05 nm, 
the barrier height has the value 0.0935 kcal/mol. The classical rate at T — 0.01 K is 1.55 x io~ 2038 sec" 1 and the 
quantum enhancemenr factor T = 1.32 x 10 9 . Even with this enhancement, the net rate is far too small to be observed. 
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Comparison between Classical and Quantum Rate at T=0.01K 



0.0 



L=L, 



-3.2508 -3.2506 -3.2504 -3.2502 -3.2500 -3.2498 -3.2496 
(L-L )/ Angstrom 



FIG. 14: Plot of logarithm of rate of crossing from one buckled state to the other using both classical TST and quantum TST, 
for a multiwalled carbon nanotube of dimensions Lq = 500 A, d\ = 10 A, d,2 = 50 A. For this nanotube the first two buckling 
instabilities occur at Lf — L = —3.2497 A and L s — Lq = —13.5555 A respectively. The first buckling instability is shown by 
an arrow and the second buckling instability is not shown in figure. The solid line is the result using classical TST, dotted line 
is the result using quantum TST. 

But following the idea of Blencowe [7j , the system which is initially in thermal equilibrium near the bottom of the 
well, can be excited to an energy level below the barrier maximum, then it will be more likely to tunnel through the 
barrier than being thermally activated over it. In such a case the rate analysis using only the reactive mode will be 
underestimating the rate by a factor of roughly 10 9 . 
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Dimensions 


Temperature 


Range of compression over which classical 
rate varies between 10 5 and 10 9 sec^ 1 


r 


Si rod 








L=1000A 








w=200 A 


300 K 


0.14 A 


» l 


d=100 A 


0.01K 


0.001 A 


« l 


Si rod 








L=20000 A 








w=200 A 


300 K 


0.6 A 


» l 


d=100 A 


0.01K 


0.006 A 


« 10 


Si rod 








L=500 A 








w=20A 


300 K 


lA 


» 1 


d=10A 


0.01 K 


0.015 A 


« 10 6 


C nanotube 








L= 5000 A 








di=50A 


300 K 


0.15 A 


« 1 


d 2 =100A 


0.01K 


o.ooi A 


« 10 


C nanotube 








L= 500 A 








di=loA 


300 K 


0.08 A 


« 1 


d 2 =50 A 


0.01K 


0.0006 A 


« 10 



TABLE I: Summary of results. 
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